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We investigate the spin and charge susceptibilities of the two-dimensional Hubbard model based 
upon the perturbative calculation in the strength of correlation U. For U comparable to a bare 
bandwidth, the charge susceptibility decreases near the half-filling as hole-doping approaches zero. 
This behavior suggesting the precursor of the Mott-Hubbard gap formation cannot be obtained 
without the vertex corrections beyond the random phase approximation. In the low-temperature 
region, the spin susceptibility deviates from the Curie- Weiss-like law and finally turns to decrease 
with the decrease of temperature. This spin-gap-like behavior is originating from the van Hove 
singularity in the density of states. 



I. INTRODUCTION 

The Hubbard modcfl has been studied extensively as 
the fundamental model for strongly correlated electron 
systems which exhibit the Mott-Hubbard metal-insulator 
transition (MIT)El. Especially, the two-dimensional Hub- 
bard model (2DHM) has been attracted much attention 
in relation to high-temperature superconductors (HTSC) 
of which normal metallic phase is quite "anomalous" in 
the sense that its property seriously deviates from that 
predicted by the conventional Fermi-liquid theory. Thus 
in this paper, we investigate 2DHM putting stress on the 
correlation effects on the charge and spin susceptibili- 
ties and grasp the signal of the "anomalous" behavior in 
the charge and spin responses near MIT from the weak- 
coupling perturbation approach. 

As for the charge degree of freedom, we are interested 
in the dependence on hole-doping 5 of the charge suscep- 
tibility. Since the system becomes incompressible near 
MIT, it is generally expected that the charge suscep- 
tibility decreases and eventually vanishes when we ap- 
proach MIT by putting 6 as zero. In fact, such a behav- 
ior has been obtained in some numerical studies on tha 
basis of the quantum Monte Carlo (QMC) methodBcl 
A quite different .result, however, has been reported in 
QMC simulationsO: The charge susceptibility diverges 
like ~ 1/8 as 6 approaches zero. The difference between 
two results is not brought about by the scatter of the 
QMC data among the groups, but is originating from 
the physical interpretation of data. Thus it is desirable 
to investigate the doping dependence of the charge sus- 
ceptibility with use of other methods. In this paper, 
by using the second-order perturbation theory (SOPT) 
with respect to the on-site Coulomb interaction U, we 
investigate the tendency of the change of the charge sus- 



ceptibility from that in the non-interacting case due to 
correlation effects. 

On the other hand, as for an anomalous behavior of the 
spin response, we would like to note the spin-gap-like be- 
havior observed in the nuclear spin relaxation rate T^ 1 of 
several HTSC materials. B~EJ In the higher-temperature 
region than the superconducting transition temperature 
T c , (TiT) -1 increases in proportional to 1/(T + 9) with 
the decrease of T, indicating that (TiT) -1 obeys the 
Curie- Weiss law. Here 9 is the Weiss temperature. As. 
T is decreased, in such materials as a Lai_ 2 ;Sr x Cu04El 
and YBa 2 Cu 3 7 _ 5 with T c ~ 90K,Q<eI (TiT)" 1 con- 
tinues to increase untill there occurs the transition to 
superconducting state. However, in such materials as 
YBa 2 Cu 3 6 . 6 with T c ~ 60K,lrE] YBa 2 Cu 4 8 ,ErEl and 
Bi 2 Sr 2 CaCu0 8 ,EJ (TiT) -1 begins to deviate from the 
Curie- Weiss law even for T > T c , and turns to decrease 
with the decrease of T. Since this decrease of (TiT) -1 
seems to be associated with the formation of the spin- 
singlet state, such a behavior has been frequently called 
"the spin gap" . Some authors have argued the origin of 
the spin-gap formation such as the spinon pairingE£l or 
bi-layer couplingjlj but it should be clarified whether it 
is possible or not to reproduce the spin-gap behavior in 
the framework of the Fermi-liquid theory without bi-layer 
coupling. From the analysis of the normal state of HTSC 
on the basis of the Fermi-liquid theory in due consider- 
ation of anti- ferromagnetic (AF) spin fluctuations ,L3 the 
following relation holds in HTSC: 



(TiT) -1 oc Xs(Q,0), 



(1.1) 



where x s (Q, 0) is the staggered spin susceptibility in the 
static limit with Q = (7r,7r). Thus in this paper, by 
calculating Xs(QiO) m SOPT, we make a qualitative in- 
vestigation on the temperature dependence of (TiT) -1 
of 2DHM. 
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By using the perturbati^calculation, several authors 
have studied 2DHM so farJiT EJ but the spin and charge 
susceptibilities which reflect the correlation effects on the 
spin and charge excitation have not been studied well 
even from the view point of SOPT. Although the pertur- 
bation method is not appropriate for the description of 
phase transition such as MIT, it is useful to study how 
correlation effects modify the behavior of these quantities 
in the Fermi-liquid phase and cause "anomalous" behav- 
iors. We can point out such a merit of the perturbative 
approach that it is free from size effects compared with 
numerical simulations and suitable for the description of 
the low-energy properties. On the other hand, there ex- 
ists such a demerit that the perturbative calculation will 
fail down for enough large value of U . However, as is 
seen in the perturbation study of the single impurity An- 
derson model, the expansion in terms of U is expected 
to converge asymptotically unless the spontaneous sym- 
metry breaking occurs, leading to that the truncation of 
the expansion in second order may give a qualitatively 
correct description of some physical properties. More- 
over, SOPT is helpful to grasp some features in the limit 
of large value of U. For example, it properly predicts 
the location of the Mott-Hubbard band at half-filling. 
In this sense, we may expect that our results on the spin 
and charge susceptibilities are qualitatively valid even for 
large value of U near half-filling. Thus it is meaningful 
to study the strong correlation effects within SOPT. 

In order to evaluate the spin and charge susceptibili- 
ties properly in the vicinity of the half-filling, we take into 
account the vertex corrections which are not included in 
the random phase approximation (RPA). These correc- 
tions are considered to be important in two-dimensional 
systems, because the charge and spin fluctuations which 
destroy long-range order are large and a simple mean- 
field description such as RPA will not be valid. Actually, 
the vertex corrections beyond RPA are important for the 
description of the charge excitations near the half-filling 
as we will see later. The calculation is performed in the 
framework of the well-known Green's function formalism. 
In order to carry out the complicated numerical computa- 
tion, we exploit the Fast-Fourier- Transformation (FFT) 
algorithm. 

The organization of this paper is as follows. In Sec. II, 
we develop the formulation to perform our calculation 
on the basis of the microscopic Fermi-liquid theory. In 
Sec. Ill, the numerical results for the charge susceptibility 
are given. We investigate the usual Hubbard model only 
with nearest-neighbor hopping terms as well as the model 
with next nearest-neighbor hopping ones which deforms 
the shape of the Fermi surface and changes the location of 
the van Hove singularity (VHS). For U comparable to the 
nearest-neighbor hopping energy, the charge susceptibil- 
ity decreases as 5 approaches zero due to the suppression 
effect of the vertex corrections. In Sec. IV, we show our 
numerical results for the dependence on S as well as T 
of the spin susceptibility. In the low-temperature region, 
the spin susceptibility shows spin-gap-like behavior due 



to VHS in the density of states (DOS). Finally in Sec. V, 
we summarize this paper and briefly comment on the re- 
lation between our results and NMR experimental results 
in HTSC. We will employ units in which both h and k& 
are taken to unity. 



II. FORMULATION 

A. Model Hamiltonian 

Let us consider tight-binding electrons on the two- 
dimensional square lattice. The non-interacting part of 
the Hamiltonian consists of the hopping term between 
nearest- and next nearest-neighbor sites. We take into 
account the on-site Coulomb repulsion U. Thus, we ob- 
tain the following model Hamiltonian: 



k<r 



(U/N)^2 E c k-qT c k<+ q i c k'l c kT, 



(2.1) 



where Cko- is the annihilation operator of an electron with 
momentum k and spin a, fj, the chemical potential, U the 
on-site Coulomb interaction and N the number of sites. 
The electronic dispersion relation is given by 



-t(cos k x + cos k y ) — t cos k x cos k y 



(2.2) 



where t and t' are, respectively, proportional to nearest 
and next-nearest neighbor hopping integrals. 



B. Self-energy in SOPT 

In general, the one-electron Green's function G(k) is 
expressed as 



G(k) = 



1 



le n + fi- E k - E(k) '' 



(2.3) 



where S(fc) is the electronic self-energy, k denotes 
(k, ie n ), e n = TtT(2n+ 1) with an integer n is the fermion 
Matsubara frequency at a temperature T, and /i is de- 
termined through the relation 



= 2^G(fc)e"" J) . 



(2.4) 



Here n is the number density of electron, r] is the positive 
infinitesimal, and we use such a shorthand notation as 



E^EE- 



(2.5) 



k n k 

With use of n, the hole-doping 5 is defined as 

5=l-7i. (2.6) 
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The Hartree term of the self-energy is expressed as 

S H = (7^G(fc)e ie "", (2.7) 

which is equal to Un/2, leading to a constant energy shift. 
Thus, we define the Green's function in the Hartree ap- 
proximation, Go(fc), as 



Go (A) 



1 



Mo - E\t ' 



(2.8) 



where fj,Q is determined through the relation Eq. (2.4) in 
which G is replaced with Go- In this number-conserving 
approximation, the Hartree term is automatically taken 
into account in Go- Therefore in the following, we do not 
include the Hartree term explicitly in each order term. 

Since the Hartree term is absorbed into /io, we need to 
calculate only the second-order term of the self-energy, 
given by 



E(fc) = U 2 Y,Xo(k~k')G (k') 

k' 

where Xo{q) is defined by 



(2.9) 



(2.10) 



Here we use another shorthand notation as q = (q, iuiv), 
where uj v = 2-kTv with an integer v is the boson Mat- 
subara frequency. 

Let us evaluate \i within SOPT, expressed as 



(2.11) 



where Sfj, is the chemical p ot enti al due to th e sec ond- 
order term. From Eqs. (EDjj), (|]|, and ( |2.1l| ), up 



to second order in U , n in Eq. ( |2.4[) is expressed as 
n = 2^2 G {k)e 1 ^ + 2^ G 2 a {k) [E(jfe) - fyt] . (2.12) 

k k 

Since the electron number is conserved, Sfi is given by 
5^ = ^r^G2(fc)S(fc). (2.13) 

In the following, we rewrite S(fc) — dfi as E(fc). 

C. Charge and spin susceptibilities in SOPT 

We define Xs(l) and Xc(q) as 



Xc(q) = [ dTe^- T (T T [ An q (r)An_ q ]), (2.14) 
Jo 



Xs(q) = (5Mb) 2 [ 1/T dTe^{T T [ S*(r)S* q ]), (2.15) 
Jo 

where and An q are, respectively, given by 

S* = (1/2) $>£ +qT Ck T - 4 +qi c ki ], (2.16) 



and An q = n q — (rt q ) with 



= £4 



k+q<T CkCT 



(2.17) 



Here g is the electron g-factor and /j,b is the Bohr mag- 
neton. The magnetic unit is fixed as <?/xb/2 = 1. 
Equations ( [2.14D and fl2.15| ) can be recast into 



Xc(q) = 2\x n (q) + X U (q)}, 



and 



Xs (q)=2\x n (q)-X n (q)], 
respectively, where x a<T (q) ls defined as 

fl/T 



(2.18) 



(2.19) 



Jo 



£4- 



qo-^ptr/ 



(2.20) 



where the subscript "conn" denotes the operation to take 
into account only con necte d diagrams. Expanding the 
right-hand side of Eq. ( 2.2C| ) with respect to U up to sec- 
ond order in the framework of the conventional perturba- 
tion theory, we obtain x aa (q) diagra mma tically shown 
in Fig. 0. Combining Eqs. (2.18) and (2.19) with the ex- 
plicit expressions for x aa (?)> we obtain Xs{q) and Xs(q) 
as 

Xc(q) = 2[xo(q)-Uxt(q) + U 2 X 3 (q) 

- 8 X (q) - U 2 {2 X corn(q) + XcoMq)}] , (2-21) 



and 



Xs(q) = 2[xo(q) + U X 2 (q) + U 2 xUq) 
- 5 X (q) + U 2 

Xcorr2 



(2.22) 



and 



where both Xcorri(<?) and x C orr2(<z) indicate the vertex 
corrections which are not included in RPA and Sx(q) is 
the correction term due to the insertion of the self-energy. 
The expressions for Xcorri(g), Xcorr2(g), and Sx(q) are 
given in the following: 

Xcorri(g) = ^2 G (fci)G (fci +q)x.o{k\ - k 2 ) 

ki,k 2 

xG Q (k 2 )Go(k 2 +q), (2.23) 
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Xcorr2 

fei ,k 2 

x G Q (k 2 )G (k 2 -q), (2.24) 
with cj>o(q) defined by 

fa(q) = -Y,Go(k)G (q-k), (2.25) 

k 

and 

Sx(q) = G o(fc)S(fc) [G (fc + q) + G (k - q)] . (2.26) 

k 

For the comparison, we also consider the susceptibilities 
calculated in RPA-like approximation up to second order 
in U, which are given by 



and 



*r A (?) = 2[xo(g) - Uxt{q) + U 2 xt{q% (2.27) 



XT A (?) = 2[xo(«) + tfxS(?) + ^X8(9)]- (2-28) 



In the following, we call this approximation "second- 
order RPA (SO-RPA)" . 



D. Method for numerical calculation 

When we carry out sum over the momentum as well 
as the Matsubara frequency in the numerical calcula- 
tion, the first Brillouin zone (FBZ) is divided into fine 
meshes and the frequency sum is terminated at a cut-off 
energy oj c . In this paper, we perform the calculation on 
a 64 x 64 lattice from the limitation of computer mem- 
ory. We should note an existence of a lower limit of tem- 
perature in which reliable results arer-abtained: As has 
been pointed out by Serene and HessE3, with a 64 x 64 
lattice, the deviation from the Fermi-liquid behavior be- 
comes significant at T/t = 0.0075 for our definition of t. 
Thus, our calculation will be carried out only for T larger 
than 0.01*. 

As for the frequency cut-off, we always choose uj c as 
64*, irrespective of T. Note that cu c is equal to 16W, 
where W is a bare bandwidth given by 4*. Such a large 
value for the cut-off energy is considered to be sufficient 
to converge the frequency sum.EHI 

We briefly explain a means to perform the compli- 
cated sum over k efficiently. The point is to arrange the 
sum into the-convolution form convenie nt for the FFT 
algorithm.EM3 For example, Xcorri in Eq. ( [2.23 ) is rewrit- 
ten as 



Yc 



rl (q) = Ci(k 1 ;q) X o(ki - fc 2 )G 1 (fc 2 ;q), (2.29) 



and 



C 1 (k- 1 q) = J2^ ky C 1 (y;q), 



respectively, Xcorri is expressed as 

Xcorri (q) oc J ^2C 1 (x;q)xo(x)C 1 (-x;q), 



(2.30) 



(2.31) 



(2.32) 



except for the normalization. Let us estimate the numer- 
ical steps necessary for the calculation with and without 
FFT. Assume that we need M steps to carry out sum 
over k. When we calculate Xcorri in Eq. ( 2.23 ) without 
FFT for a fixed q, the steps of the order of M 2 are neces- 
sary. On the other hand, if we exploi t the FFT algorithm 
for the calculation of Xcorri in Eq. (2.32), we need only 
the steps of the order of Mlog 2 M to finish all the sums. 
This difference leads us a drastic reduction of the task in 
the numerical calculation. 



III. CALCULATED RESULTS FOR CHARGE 
SUSCEPTIBILITY 

A. Doping dependence 

Now we investigate the dependence on 5 of the uniform 
charge susceptibility in the static limit, Xc- We first con- 
sider the case only with nearest-ne ighb or hopping in the 
dispersion relation given by Eq. (2^2). The results for 
this case are shown in Fig. p|. For U = 0, Xc increases 
logarithmically due to VHS, as 6 approaches zero. In 
the interacting case, \ c is totally suppressed both in SO- 
RPA and SOPT. For U comparable to t, Xc still increases 
with the decrease of 5 in SO-RPA, but in contrast, it 
turns to decrease near the half-filling in SOPT as 5 ap- 
proaches zero. The large suppression of Xc in SOPT near 
the half-filling is due to the vertex corrections which are 
not included in RPA. 

In order to discuss the physical meaning of our calcu- 
lated results, we would like to comment on the relation 
between the formula for Xc of the Landau's Fermi-liquid 
theory and one of the microscopic Fermi-liquid theory 
developed by LuttingerBj In the Landau's Fermi-liquid 
theory, Xc is phcnomcnologically expressed as 



Xc 



(3.1) 



where C\{k; q) = Go(k)Go(k + q). When we symbolically 
represent the Fourier transform of Xo(q) and C\{k; q) as 



where m* is the renormalized mass, Xc the bare suscep- 
tibility, and Fq the Landau's parameter. Based upon the 
Luttinger's formulation, we obtain another expression for 

Xc as 
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! _ dRe£(k,0) 



(3.2) 



B. Effects of next- nearest neighbor hopping 



where is the dispersion relation for the quasi-particle 
given by 

££ = z k [£ k -M + Re£(k,0)], (3.3) 
and Zk is the renormalization factor determined by 
<9Re£(k,w) 



^ = 1- 



dui 



(3.4) 



Comparing Eq. ( ft. 2] ) with Eq. ( |3.l| ), we notice approxi- 
mately the following correspondence: 



and 



<9Re£(k,0) 
9/i 



(3.5) 



(3.6) 



From Eq. (3.2), we see that Xc is nothing but the aver- 
age over the Fermi surface of the factor 1 — <9ReI]/<9/x. 
It does not depend exp licitly on the renormalized mass 
in contrast to Eq. (3.1), but is determined by the bare 
DOS and the factor 1 — <9Re£/<9/x which is incorporated 
into our calculation as the vertex corrections. Thus the 
correlation effect on Xc is essentially due to the factor 
1 — dKeH/d/i. Note that Xc is not enhanced by the renor- 
malized mass m* as naively expected from Eq. (3T). This 
is contrasted with the one-dimensional Luttinger liquid 
where the charge susceptibility is determined solely by 
the renormalized DOS.E-3 

Our numerical results imply that 1 — <9Re£/<9/i tends to 
be zero as /x approaches the half-filling value, suggesting 
that Xc vanishes continuously if we can describe the MIT 
transition properly. Our results are consistent with that 
of the QMC simulations by Moreo et. alu and Oteukao 
but not with the assertion by Furukawa and ImadaQ that 
the charge susceptibility diverges in the vicinity of the 
half-filling. In the framework of the Fermi-liquid theory, 
the increase of the charge susceptibility due to correla- 
tion effects does not occur as discussed above. Thus the 
conclusion in Ref. 5 implies that some non-Fermi-liquid 
state realizes in the metallic phase near MIT. However, 
the results in the infinite-dimensional Hubbard model in- 
dicate that there is no intermediate exotic state between 
the FermijJjquid metallic phase and the Mott insulat- 
ing phase.Ej Since no symmetry breaking occurs at finite 
temperature in two-dimensional systems, we expect that 
there is no phase transition from the Fermi-liquid state 
to a non-Fermi-liquid one as n approaches unity. Thus 
our results based upon the Fermi-liquid theory plausi- 
bly may describe the qualitative behavior of the charge 
susceptibility near MIT. 



Since Xc depends on the bare DOS, it is sensitively af- 
fected by the band structure. In this subsection, we con- 
sider the effect of the next-nearest hopping term which 
changes the location of VHS and deforms the shape of 
the Fermi surface. The results for the case of t'/t = 0.2 
is shown in Fig. [}| It can be seen that Xc increases mono- 
tonically as S approaches zero, indicating that we cannot 
obtain the suppression effect large enough to overcome 
the enhancement of the bare DOS due to VHS. In this 
case, VHS is located in the region of oj > 0, where uj is 
the energy measured from the Fermi level. Thus the bare 
DOS at the Fermi level is always too small to produce 
strong correlation effects through the vertex corrections 
when we decrease hole-doping by changing the chemical 
potential. The suppression effect due to the vertex cor- 
rections is weak when the bare DOS at the Fermi level is 
not large enough. This tendency can be confirmed from 
the results for the case of t'/t = —0.2, which are shown 
in Fig. |J. In this case, VHS is located in the region of 
uj < 0. When S is decreased with the change of fi, the 
Fermi level inevitably crosses the position of VHS at some 
value of \x. At the corresponding hole doping, Xc shows 
a shallow minimum for U/t — 1.0. This minimum is due 
to the maximum of the suppression effect of the vertex 
corrections. Thus we conclude that the large bare DOS 
at the Fermi level results in the strong suppression effect 
of Xc through the vertex corrections which overcomes the 
enhancement effect of Xc due to the bare DOS, leading 
to the total decrease of Xc in the vicinity of MIT. 



IV. CALCULATED RESULTS FOR SPIN 
SUSCEPTIBILITY 



A. Doping Dependence 



Let us discuss the numerical results for the static spin 
susceptibility at q = Q, Xs- The dependence on S of Xs 
is depicted in Fig. 0. Because of the nesting properties of 
the Fermi surface, the AF correlation develops near the 
half- filling and results in the enhancement of Xs- We also 
show the results for the case in the presence of the next 
nearest-neighbor hopping t' in Figs. |^ and [?]. In both 
cases of t' — 0.2 and t' = —0.2, the nesting properties 
of the Fermi surface is much weakened. Thus the en- 
hancement of AF correlation is reduced compared with 
the case for t' = 0, but there is no qualitative difference 
in the doping dependence between these cases. 
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B. Temperature Dependence: Spin-gap-like Behavior 

In this subsection, we investigate the temperature de- 
pendence of Xs- The results for the case with only 
the nearest-neighbor hopping are shown in Fig. ||. The 
electron number is chosen as n — 0.94. In the high- 
temperature region, \s follows the Curie- Weiss-like be- 
havior. On the other hand, in the low-temperature re- 
gion, the spin-gap- like behavior manifests. Since it is 
observed even in the case of U = 0, this spin-gap-like 
behavior is essentially due to the band effect: In the 
non-half-filled case, VHS's at k = (±7r, ±7r) are not on 
the Fermi surface and then the contributions from VHS 
to the spin susceptibility decrease as the temperature is 
lowered. Consequently, \s shows the spin-gap-like be- 
havior in the low-temperature region. Although the spin- 
gap-like behavior is very small for U = near the half- 
filling, it is much enhanced due to the correlation effect 
and results in the large spin-gap-like behavior as shown 
in Fig. ||. Note that such enhancement is not brought 
about by the speciality of RPA, because we also obtain 
the robust spin-gap-like structure even with the vertex 
corrections beyond RPA. 

When we introduce the next nearest-neighbor hopping, 
the spin-gap-like behavior disappears as shown in Figs. O 
and [l0| for t' = 0.2 and t' = —0.2, respectively, with 
n = 0.94. The next nearest-neighbor hopping changes 
band structure and the effect of VHS becomes too small 
to give rise the spin-gap behavior in these cases. For 
the case of t' = —0.2, the large enhancement in the low- 
temperature region can be understood as the effect of 
VHS, because the location of VHS is very close to the 
Fermi surface. At the lower temperature than O.Olt, it is 
expected that the spin-gap behavior may manifest, but 
we do not carry out the calculation in this temperature 
region because of the lack of the reliability of our cal- 
culation on a 64 x 64 lattice as well as the limitation 
of computer memory. Instead of calculating \ s at lower 
temperature, we obtain the temperature dependence of 
Xs at n = 0.86. The results are shown in Fig. If. Since 
the difference between the position of VHS and that of 
the Fermi level is O.Obt at this filling, the spin-gap-like be- 
havior is observed even in the present temperature range 
in contrast to the case of n = 0.94. Thus our spin-gap- 
behavior is originating from VHS, not the nesting prop- 
erty of the Fermi surface. 



V. SUMMARY AND DISCUSSION 

We have investigated the doping and temperature de- 
pendence of the charge and spin susceptibilities by using 
SOPT. The charge susceptibility is much suppressed due 
to the correlation effect through the vertex corrections 
in the vicinity of the half-filling for U comparable to t. 
It has been explicitly shown that the vertex corrections 
beyond RPA are indispensable in order to obtain this be- 



havior. The suppression of the charge susceptibility near 
the half-filling indicates a precursor of the formation of 
the Mott-Hubbard gap. Such a behavior of the charge 
susceptibility will be qualitatively correct as long as the 
system is kept in the Fermi-liquid state, but we cannot 
perfectly deny the possibility that the charge susceptibil- 
ity jumps abruptly from a finite value to zero or diverges 
at the half-filling if the system turns to be in a non- Fermi- 
liquid state due to the speciality of two-dimensional sys- 
tems in combination with the strong correlation effect. In 
this sense, we need further studies beyond perturbative 
calculation in order to understand the behavior of the 
charge susceptibility near MIT in the two-dimensional 
Hubbard model. 

In the study of the spin susceptibility, we have obtained 
the spin-gap-like behavior not only the case for t' = but 
also the case with non-zero next nearest-neighbor hop- 
ping. The origin of our spin-gap-like behavior is ascribed 
to the band effect such as VHS in combination with the 
enhancement effect due to electron correlation: In the 
static spin susceptibility at q = (ir, n), VHS gives rise to 
a small spin-gap-like structure, which will be much en- 
hanced by correlation effects. We have emphasized that 
such an enhancement is not brought about by the special- 
ity of RPA. Even if we take into account the diagrams for 
the vertex corrections beyond RPA, the same structure 
has been obtained. Those vertex corrections are impor- 
tant in the two-dimensional model in which quantum and 
thermal fluctuations destroy any long range order. 

Finally, we give a brief comment on the spin gap in 
NMR experiments on HTSC.ETL3 It has been argued that 
the origin of the spin-gap-like behavior observed in those 
experiments is the condensation of paired spinons in the 
RVB pictureEj or the bi-layer coupling in the Fermi-liquid 
picture£j As we have seen in this paper, even in the 
Fermi-liquid theory without bi-layer coupling, the spin- 
gap-like behavior has been reproduced by the effect of 
VHS combined with the enhancement effect due to elec- 
tron correlation. Although it is difficult to adjust quan- 
titively the size of the gap and the band structure in our 
calculation to the experimental results, there are some 
indirect evidences which imply the relevance of VHS to 
the spin-gap behavior. In the NMR experimental result 
on YB2Cu307_iE3, it can be seen that the spin-gap be- 
havior becomes weak as S approaches zero. It has been 
pointed out that the DOS of YB 2 Cu 3 07_5 has a peak 
near the Fermi level due to VHS for i5 = and-jthis peak 
runs away from the Fermi level as 5 increasesca. Thus it 
may be plausible to attribute the spin-gap-like behavior 
of this material to VHS in DOS enhanced due to strong 
electron correlation. 
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FIG. 1. Diagrams for x"" within second order with re- 
spect to the on-site Coulomb interaction U. The solid and 
the dashed lines, respectively, indicate the bare Green's func- 
tion Go and U. Note that the Hartree term is included in the 
Go-line. 



FIG. 2. Charge susceptibility plotted as a function of the 
hole-doping 8 for t' = and T — 0.04i. Open squares linked 
by thin solid line indicates the results for U = 0. The results 
for SO-RPA and SOPT for U — t are denoted by open circles 
connected by broken line and solid circles linked by thick line, 
respectively. 



FIG. 3. Charge susceptibility plotted as a function of the 
hole-doping S for t' = 0.2* and T = 0.04*. The results are 
shown in the same lines and symbols as in Fig. 2. 



FIG. 4. Charge susceptibility plotted as a function of the 
hole-doping S for *' = -0.2* and T = 0.04*. The results are 
shown in the same lines and symbols as in Fig. 2. 
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FIG. 5. Spin susceptibility as a function of the hole-doping 
S for t' = and T = 0.04*. Open squares linked by this 
solid line indicates the results for U — 0. The results for 
SO-RPA and SOPT for U = 0.5* are denoted by open circles 
connected by broken line and solid circles linked by thick line, 
respectively. 

FIG. 6. Spin susceptibility as a function of the hole-doping 
S for *' — 0.2* and T = 0.04*. The results are shown in the 
same lines and symbols as in Fig. 5. 

FIG. 7. Spin susceptibility as a function of the hole-doping 
S for *' = -0.2t and T = 0.04*. The results are shown in the 
same lines and symbols as in Fig. 5. 

FIG. 8. Spin susceptibility plotted as a function of the tem- 
perature T for t' = and n = 0.94. Open squares linked by 
this solid line indicates the results for U — 0. The results for 
SO-RPA and SOPT for U = 0.5* are denoted by open circles 
connected by broken line and solid circles linked by thick line, 
respectively. 

FIG. 9. Spin susceptibility plotted as a function of the tem- 
perature T for t' = 0.2* and n = 0.94. The results are shown 
in the same lines and symbols as in Fig. 8. 

FIG. 10. Spin susceptibility plotted as a function of the 
temperature T for *' = —0.2* and n = 0.94. The results are 
shown in the same lines and symbols as in Fig. 8. 

FIG. 11. Spin susceptibility plotted as a function of the 
temperature T for *' = —0.2* and n = 0.86. The results are 
shown in the same lines and symbols as in Fig. 8. 
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